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The production of cc and bb is studied in pp collisions at = 7 TeV with the LHCb detector. The results of 
these measurements are compared with different theoretical models. Results and prospects are also shown for 
exotics: the mass measurement of the X(3872) and the search of the X(4140). Most of the presented results 
are based on the 2010 LHCb dataset (35 pb _1 ). 



1. Introduction 

Heavy quarkonium production remains a challenging problem for the understanding of Quantum Chromo- 
dynamics (QCD). At the centre-of mass energies for proton-proton collisions at the Large Hadron Collider, cc 
pairs are expected to be produced predominantly via Leading Order gluon gluon interactions, which can be 
computed using perturbative QCD, followed by the formation of the bound charmonium states described by 
non perturbative models. Recent approaches make use of non-relativistic QCD factorization (NRQCD) which 
assumes a combination of the colour-singlet and colour-octet cc as it evolves towards the final bound state via 
the exchange of soft gluons [1]. Next-to-Leading Order (NLO) QCD corrections in charmonium and bottomo- 
nium production are also essential for the description of the experimental data [2, 3]. Studies of i/j(2S), T(1S) 
and Xc production cross sections are presented. 

The X(3872) meson is an exotic meson discovered in 2003 by the Belle collaboration in the — X(3872)K ± , 
X(3872) — >• J/^7r + 7r _ decay chain Q. Its existence was confirmed by the CDF 0, DO @ and BaBar Q 
collaborations. The X(3872) particle, together with several other new states subsequently observed in the mass 
range 3.9-4.7 GeV/c 2 , has led to a resurgence of interest in exotic meson spectroscopy [8]. Several properties of 
the X(3872) have been measured; however, its nature is still uncertain and several models have been proposed. 
First, it is not excluded that the X(3872) is a conventional charmonium state with one candidate being the 
T] C 2(1D) meson Q. However, the mass of this state is predicted to be far below the observed X(3872) mass. 
Given the proximity of the X(3872) mass to the D*°D° threshold, one possibility is that the X(3872) is a 
loosely bound deuteron-like D*°D° molecule, i.e. a ((uc) - (cu)) system [8|. Another more exotic possibility is 
that the X(3872) is a tetraquark state [9[. 

The CDF experiment has reported > 5 cr evidence for X(4140) state (also referred to as y(4140) in 
the literature) in 6.0 fb _1 of pp data collected at Tevatron [10]. The relative rate was measured to be 
B(B+ X(4U0)K+) x £(X(4140) J/V^)/£(£ + J/^<j>K + ) = 0.149 ± 0.039 ± 0.024. This observation 
has triggered wide interest among model builders of exotic hadronic states. It has been suggested that the 
X (4140) could be a molecular state [U, G3, a tetraquark state [HI, [H|, a hybrid state [HI, [l6| or a re-scattering 
effect E 03. 

2. \I>{2S) Production Cross Section 

Two decay modes of the ip(2S) meson have been studied [l9j: ip(2S) —> /i + ji~ and ip(2S) —> J/'0(/i/i)7r + 7r _ . 
The differential cross-section for the inclusive i/j(2S) production is computed as 

dp T KPT! L int e(p T )B^(2S) f)A PT {> 

where A^( 2 5)(pt) is the number of observed ip(2S) decays, L int is the integrated luminosity, c(jpt) is the total 
detection efficiency including acceptance effects, B(ip(2S) — )■ /) is B(ip(2S) — » J/i/j7r+7r~)B(J/'ip — >> for 
ip(2S) — >• J/i/j(jHjii)7T^7T~ and the dielectron branching ratio (B(ip(2S) — >• e + e - )) for ip(2S) — >> and Apr = 

1 GeV/c is the bin size. The dielectron branching ratio is used, assuming lepton universality, since it has a much 
smaller error than the dimuon one (2.2% and 10% respectively). In order to estimate the number of ip(2S) signal 
events, a fit is performed independently in each (p?) bin. In Fig. [1] a comparison between the measurements 
presented of the differential ^(2S) production cross-section multiplied by the ip(2S) —> /i + \i~ branching ratio 
and a recent theory prediction [20j for prompt ip(2S) in the rapidity range 2 < y = \lTi ^[_ p p z < 4.5 is shown. 
Here E and p z are the ip(2S) energy and momentum in the z direction measured in the pp centre-of-mass frame; 
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Figure 1: Comparison of the LHCb results for the differential production cross-section of ij)(2S) with the predictions for 
prompt production by a NLO NRQCD model [20||. LHCb data include also ip(2S) from b. 



the z-axis is defined along the beam axis in the LHCb frame, oriented from the VELO to the Muon Detector. 
The differential cross-section in [20| has been evaluated for the i/j(2S) prompt production at the LHC at next- 
to-leading order in non-relativistic QCD, including color-singlet and color-octet contributions. The integrated 
cross-section in the full range of pr and y respectively to ijj(2S) —> fi + fi~ and i/j(2S) —> J/^7r + 7r~ is found to 
be: 

cr(0 < p T < 12 GeV/c, 2 < y < 4.5) = 1.88 ± 0.02 ± 0.31±g;|| fib 

cr(3 < p T < 16 GeV/c, 2 < y < 4.5) = 0.62 ± 0.04 ± 0.12±g;^ fib (2) 

where the first uncertainty is statistical, the second is systematic and the third is the uncertainty due to the 
unknown polarization. 



3. T(1S) Production Cross Section 

T(16') meson has been studied in the decay mode T(1S) — » fi + fi~ [2l[. The double differential cross-section 
for the inclusive T(15') production is computed as 

d 2 o~ , ) = N r{ls) (p T ,y) 

dp T dy [PT ' y) L int e(p T ,y)B(T(lS) -+ fi+fi-)Ap T Ay [ ] 

where Nx(is)(PT,y) is the number of observed T(1S) — >• fi + fi~ decays, e{pT-,y) is the total detection efficiency 
including acceptance effects, Li nt is the integrated luminosity, B(T(1S) — > /i + fi~) is the branching fraction, 
and AprAy = 1 x 0.5 GeV/c is the rapidity and pr bin sizes. In order to estimate the number of Y(15) signal 
events, a fit is performed independently in each of the 15 pr times 5 y bins. The double differential cross-section 
as a function of pr and y is shown in Fig. [2l The integrated cross-section in the full range of y and pr is found 
to be 

a(pp -> T(1S)X; p T (T(lS)) < 15 GeV/c; 2 < y(?(lS)) < 4.5) = 108.3 ± 0.7l^;g nb (4) 

where the first uncertainty is statistical, and the second systematic. The latter includes ~t\% 8 nb from the 
unknown polarization, ±10.8 nb from the luminosity determination and ±22.0 nb from other sources. The 
integrated cross-section is about a factor 100 smaller than the integrated J/ip cross-section in the identical y 
and pt region [22], and a factor three smaller than the integrated T(15) cross-section in the central region \y\ < 2 
as measured by CMS [23]. The CMS and the LHCb measurements agree; the difference in the integrated cross 
section results is due to the different rapidity ranges. Fig. [3] shows the comparison of the LHCb cross-section 
measurement in bins of y integrated over pt with the same measurement of CMS, in the y ranges covered by 
the two experiments. 
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Figure 2: Differential T(l*S f ) production cross-section as a function of pr in bins of rapidity. 
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Figure 3: Differential T(1S) production cross-section as a function of y integrated over pr as measured by the CMS [22 
and LHCb experiments. 



4- o-(xc2)/o-(xd) Production Cross Sections Ratio 

Xci and Xc2 have been reconstructed in the final state J/tpj. The production cross-section ratio of the Xc2 
and Xci states is measured Q using 



cr( X c2) _ N Xc2 g/^Tfg g( Xcl -» J/^ 7 ) 



(5) 



where B(xd ~^ J/^l) an d B(xc2 J/^l) are the Xci an d Xc2 branching ratios to the final state J/^HS ^j/^ 
(e* c2 ) is the efficiency to reconstruct and select a J/ip (7) from Xc decay and e*^ is the efficiency to select 
the x c candidate. The measurement method consists of extracting the two N Xci yields (for Xci an d Xc2) from 
an unbinned maximum likelihood fit to AM = M(xc) ~ M(J/rf) mass difference distribution. With the ratio 
of the resolution parameters and the mass differences fixed, a fit is then performed to the data in the full 
J/ip pt range J/ip pr G [3,15] GeV/c in order to extract the resolution scale cr res (Xci)- Here, the sample is 
subdivided into candidates with converted (after the magnet) and non-converted photons, in order to account 
for the different calorimeter resolution in the two cases. Converted photons are identified by the presence of 
activity in the scintillator pad detector at the entrance to the calorimeter system. Photons that convert before 
the magnet have a low probability to be reconstructed, because either one or both electrons are swept out of 
the detector acceptance by the magnetic field, and are not considered here. The resolution scales, cr res (Xci) 7 
are measured to be 22.8 ±1.1 MeV/c 2 and 18.4 ± 0.4 MeV/c 2 for converted and non converted candidates, 
respectively. The fit is then performed in bins of J/ip pr- For each bin the value of cr res (Xci) is fixed to the 
value extracted from the fit to the full range J/ip pr- The presence of Xc polarized states would modify the 
efficiencies calculated from the Monte Carlo. In order to take into account possible polarization scenarios, the 
relevant combination of weights has to be taken into account. The results from the not-converted and converted 
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Figure 4: The ratio cr(xc2)/o"(Xci) in bins of in the range of J/tp pr ' [3; 15] GeV/c. The internal error bars correspond 
to the statistical error on the % c i and %c2 yields; the external error bars include the contribution from the systematic 
uncertainties (apart from the polarization). The shaded area around the data points (black) shows the maximum effect 
of the unknown \c polarization on the result. The two other bands correspond to the ChiGen MC generator theoretical 
prediction [2^ (in blue) and NLO NRQCD [H] (in red). 



samples are combined by sampling the combination of the statistical and uncorrelated systematic uncertainties 
using a toy Monte Carlo. Finally, the correlated systematic uncertainty from the branching ratio is calculated 
using the combined central values for (t(Xc2)/&(Xci)- The preliminary result for the ratio of the prompt Xc2 to 
Xd production cross-sections as a function of J/ip pt is given in Fig. HJ Comparisons to the theory predictions 
from the ChiGen MC generator [25[ and from the NLO NRQCD calculations [26], in the rapidity range [2; 4.5], 
are also shown in the figure. Fig. [4] also shows the maximum effect of the unknown Xc polarization on the result, 
shown as a black shaded area around the data points. The upper limit of the shaded area corresponds to the 
spin state ( Xci : m J = 1; Xc2 - mj = 2) and the lower limit corresponds to the spin state ( xa • raj = 0; Xc2 
: mj = 0). The results are broadly in agreement at high J /if; pr with the colour singlet model, however, they 
are not yet precise enough to rule out sizeable colour octet terms. There are indications of a discrepancy in the 
mid to low J/tp pt region. This may be explained by a more complete modeling of the transition from the high 
to low J/ip pt regions, which is sensitive to non-pert urbative effects and/or sizeable higher-order perturbative 
corrections. 



5. X(3872) Mass Measurement 

Inclusively produced X(3872) mesons are reconstructed and selected in the X(3872) — ^ J/?/;7r + 7r _ , J/ip — » 
decay mode [27j. A momentum scale calibration is performed to account for a mixture of effects related 
to imperfections in the knowledge of the magnetic field map and of the alignment of the tracking system. 

5.1. Momentum Scale Calibration 

The momentum scale is calibrated using a large sample of J/ip — >• /i + \i~ decays. This calibration, gives an 
overall scale factor that is applied to all measurements of particle momenta. The calibration is checked to be 
also valid for two-body decays of the T(15), D° and K°. In addition, its effect on tp(2S) —> J/^7r + 7r _ , which 
have kinematics similar to the X(3872) — >> J/^7r + 7r _ , is studied. The uncalibrated value of the ip(2S) mass in 
data is 3685.94 ± 0.06(stat) MeV/c 2 . After the momentum scale calibration it becomes 3686.12 ± 0.06(stat) 
MeV/c 2 , in good agreement with the PDG value of 3686.09 ± 0.04 MeV/c 2 . The measured J/ip mass after 
alignment and calibration is checked to be stable over the whole 2010 data-taking period. 
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Figure 5: Invariant mass distribution of J/^n + n~ (black points with statistical error bars) and same-sign J/ip^^ 
(blue filled histogram) candidates. The red curve is the result of the fit described in the text. The insert shows a zoom 
of the region around the X(3872) mass. 



5.2. Results 

The masses of the ijj(2S) and X(3872) are determined from an extended unbinned maximum likelihood fit of 
the reconstructed J/^7r + 7r~ mass in the interval 3.6 < M J /^ 7r + 7r - < 3.95 GeV/c 2 . The functional form of the 
background is studied using the same-sign pion events. The ift(2S) and JT(3872) signals are each described with 
a Voigt function denned as the convolution of a non-relativistic Breit-Wigner with a Gaussian function. The 
intrinsic width of the ^(25) is fixed to the PDG value, T^ 2 s) = 0.317 MeV/c 2 . The X(3872) intrinsic width 
is poorly known. The BaBar [28] and Belle [4J collaborations have published 90% confidence level limits of 
Tx(3872) < 3.3 MeV/c 2 and 1^(3872) < 2.3 MeV/c 2 , respectively. The results of the Belle and BaBar analyses 
have been combined by CDF to yield 1^(3872) = 1.3 ±0.6 MeV/c 2 [29]. For these studies, the procedure adopted 
in [29| has been adopted and the natural width has been fixed to this value in the default fit. Fig. [5] shows the 
invariant J/^7r + 7r~ mass distributions for opposite-sign (black points) and same-sign (blue filled histogram) 
candidates. Clear signals for both the i/)(2S) and the X(3872) can be seen. From a comparison of fits with 
and without the X(3872) component the statistical significance of the X(3872) signal is estimated to be 9 a. 
The uncertainties on the parameters reported by the fit are in good agreement with expectations based on toy 
Monte Carlo studies. The preliminary LHCb result 

M x(38 72) = 3871.96 ± 0.46 ± 0.10 MeV/c 2 (6) 

is in good agreement with the published measurements 0, IU HH Hlj and with their average, 3871.56 ± 0.22 
MeV/c 2 . The new world average value including LHCb measurement, 3871.63 ±0.20 MeV/c 2 , is consistent 
within uncertainties with the sum of the D° and D*° masses, 3871.79 ± 0.29 MeV/c 2 , computed from the 
results of the global PDG fit of the charm meson masses. 



6. X(3872) Production Cross Section 

The product of the inclusive production cross-section cr(pp — >• X(3872) ± ...) by the X(3872) — » J/^7r + 7r _ 
branching fraction is computed as (30| 

jycorr 

a(pp X(3872) ± ...) x 5(X(3872) J/^tt') = ^g^l ( 7 ) 

where ^V"x(3872) ^ s ^ ne efficiency-corrected yield of X(3872) — >• J/^(/i + /i _ )7r + 7r _ signal decays, rjtot is a multi- 
plicative factor to the efficiency that accounts for known differences between the data and simulation, L int is the 
integrated luminosity, and the notation B(...) is used for branching fractions. The X(3872) signal yield is deter- 
mined from an extended, unbinned maximum likelihood fit of the reconstructed J/^7r + 7r _ mass in the interval 
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3.82 < M J ^ 7T + 7T - < 3.95 GeV/c 2 . The fitting function is the same used for the X(3872) mass measurement. 
The mass resolution is fixed in this analysis to 

a MC 

<r(X(3872)) = (7(^(25)) X{ ^J 2) = 3.26 ± 0.10 MeV/c 2 (8) 

where a(ip(2S)) = 2.48 ± 0.08 MeV/c 2 is the mass resolution fitted in data for the kinematically similar decay 
ip(2S) —> J/^7r + 7r _ , and where cr^^ 872 ^ and o-^f^s) are ^ ne mass resolutions determined from fully simulated 
Monte Carlo events. As in section [5^ the natural width is fixed to 1.3 MeV/c 2 . The effect of fixing parameters 
that are not perfectly known is investigated as part of the systematic studies. The X(3872) signal yield is 
approximately 68% of that reported in [27] due to the additional trigger requirements and fiducial cuts on y and 
Pt- Applying the procedure discussed above, and using the pt~ and y— dependent efficiency from simulation, 
the efficiency-corrected signal yield is ^x°(3872) = 9597 ±2217, where the quoted uncertainty is statistical. Two 
factors enter into the determination of the correction factor rjtot- The first is a factor 1.024 ± 0.011 which 
accounts for differences in the efficiency of the muon identification in the data and simulation. The second is a 
factor of 0.96 ± 0.02 that accounts for observed differences in the efficiency of global event cuts applied in the 
trigger. Multiplying these values together gives r] to t = 0.983 ± 0.023. The uncertainty on this number is taken 
into account in the estimation of the systematic uncertainty. The preliminary measurement is 

a(pp X(3872) + ...) x £(X(3872) J/^tt+tt-) = 4.74 ± 1.10 ± 1.01 nb (9) 

where a(pp — ^ X(3872) + ...) is the cross-section for producing an X(3872) particle in pp collisions at y/s = 7 TeV 
(either promptly or from the decay of other particles) with a transverse momentum between 5 and 20 GeV/c 
and a rapidity between 2.5 and 4.5. 



7. Search of the X(4140) 

The X(4140) has been studied using a data sample of approximately 0.376 fb _1 [33| selecting B + — » J/i/j<fiK+, 
(j) K + K~ . No narrow structure has been seen near the threshold in M{J — M(J/ip) distribution as 
shown in Fig. [6] In the CDF analysis, they fit their data with spin relativistic Breit-Wigner function on top of 
three-body phase-space, all smeared with the detector resolution [To|- To quantify the disagreement with CDF 
the same function has been used. The efficiency dependence is extracted from the MC simulations and applied 
as a correction to the three-body phase-space function. Mass and width of X(4140) peak are fixed to the values 
obtained by the CDF collaboration. The mass difference resolution is determined from the B + — >• X(4140)if + 
MC simulation. The fit to M{J/^(j>) - M(J/^) distribution gives a X(4140) amplitude of 6.9 ± 4.7 events 
(Fig. [6ti) and has a confidence level (CL) of 3%. When using a quadratic polynomial instead of three-body 
phase-space function for the background, the preferred value of the X(4140) amplitude is zero (it is restricted 
not to go below zero), with a positive error of 3 events. This fit is shown in Fig. [6b and has a confidence level 
of 11%. Using the B+ J/^(j)K + yield multiplied by this efficiency ratio (B+ X(4140)if+, X(4140)_-> 
J/^<t>)/{B + -> J/^<t>K + ) and multiplied by the CDF value for B(B+ -> X(4140)if +)/£(£+ -> J/^<t>K + ) P3, 
leads to a prediction of observed 35 ± 9 ± 6 events, where the first uncertainty is statistical from the CDF data 
and the second is systematic including both CDF and LHCb contributions. The central value of this estimate 
is illustrated in Fig. [6] The CDF result disagrees by three standard deviations with the fit to LHCb data using 
the polynomial background. The disagreement is 2.4 a when the efficiency-corrected three-body phase-space 
background shape is used. 



8. Exclusive Dimuon Production 

Exclusive particle production in proton-proton collisions are elastic processes in which the protons remain 
intact, p + p — >• p + X + p, and the additional particles are created through photon and/or gluon propagators. In 
the case of purely photon propagators, this is an electromagnetic process which can be theoretically calculated 
with high accuracy. When gluons are involved, these processes constitute an important testing ground for 
QCD, since the object that couples to the proton must be colourless. Thus the pomeron (two gluon states) 
or possibly an odderon (three gluons), predicted in QCD but never unambiguously observed, can be studied 
in a clean experimental environment. The cleanest experimental manifestation of these occur in final states 
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Figure 6: Fit of X(4140) signal on top of smooth background to M(J/ip<j)) - M(J/i/>) for the B + J/t/j^K^ data. The 
solid line represents the result of the fit to LHCb data. The dashed line on top illustrates the expected signal amplitude 
from the CDF results as explained in the text. The top and bottom plots differ by the type of the background function 
used in the fit: a) efficiency-shaped three-body phase-spaced; b) quadratic polynomial. 



containing two protons and two muons or in events containing two protons, two muons and a photon. The 
former can be produced in the diphoton process giving a continuous dimuon invariant mass spectrum, or in 
the photon-pomeron process which can produce (/>, J/i/j, \/j(2S), Y(15), T(25) or T(35) which decay to two 
muons. The latter is a signal for double pomeron exchange which produces Xc(Xb) that decay to J/ip (T) plus a 
photon. The final state protons are only marginally deflected, go down the beam-pipe, and remain undetected. 
The experimental signal therefore in LHCb is a completely empty event except for two muons and possibly a 
photon. However, because LHCb is not hermetic, there will be sizeable backgrounds from non-elastic processes 
where the other particles travel outside the detector acceptance. 

The cross-section cr, is calculated [34[ from the number, TV, of selected events having corrected for efficiency, 
e, and purity, p, and dividing by the luminosity, L, via a = (pN)/(eL). The efficiency for selecting the events 
has been determined from simulation. The measured cross-sections are: 

^j/^+m( 2 <^+'W< 4 - 5 ) = 474±12±51±92pb 

^(25)^+m( 2 <^+'W< 4 - 5 ) = 12.2 ± 1.8 ± 1.3 ± 2.4 pb 

^x c o^J/^/i+M7( 2< ^+'W^7< 4 - 5 ) = 9.3±2.2±3.5±1.8pb (10) 

cr Xcl ^ J /^ 7 ^ /i + M7 (2 < 7^+, 77/i,7? 7 < 4.5) = 16.4 ± 5.3 ± 5.8 ± 3.2 pb 

^^J/^mW 2 < V^W' 7 ^ < 4.5) = 28.0 ± 5.4 ± 9.7 ± 5.4 pb 

^^+M-p( 2< ^+'W< 4 - 5 ' M /i+/i- >2.5GeV/c 2 ) = 67±10±7±15pb 

where the first uncertainty is statistical, the second is systematic, and the third comes from the estimate of the 
luminosity. Note that these numbers are cross-section times the branching ratio into the final state of interest, 
and all final state particles are required to be between pseudorapidities of 2 and 4.5. The cross sections are 
quoted in a limited pseudorapidity range as the models under consideration have pseudorapidity dependence. 
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9. Summary 

LHCb performed many analysis of the quarkonium states using 2010 collected data (35 ptr 1 ). The measure- 
ment of the production cross sections of charmonium and bottomonium states (-0(25), T(15), Xc) are useful 
to test theoretical models. For the exotic states, the measurement of the X(3872) mass and cross section has 
been performed and the CDF narrow X(4140) has been studied using a dataset of 376 pb _1 and its existence 
not confirmed. LHCb has a very high J/ip statistic in 2011 data that will allow to have a lot of new results in 
the near future. 
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